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illustration  is  for  steady  high  Magnetic  Reynolds  Number  channel  flows. 
Here  we  exhibit  the  nonuniform  manner  in  which  the  external  magnetic 
field  is  diffused  into  the  plasma  and  the  corresponding  nonuniform 
deceleration  of  the  flow. 
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INTRODUCTION 


Advances  in  the  performance  levels  of  devices  based  upon  the 
interaction  of  electromagnetic  forces  with  a conducting  fluid  depend  in  no 
small  measure  upon  the  understanding  and  control  of  electrical  and  fluid 
nonuniformities  present  in  such  machines.  Such  nonuniformities  may  be 
spatially  induced  by  the  presence  of  boundary  walls  or  externally  imposed 
fields  (boundary  layers,  separation/inter  action  zones,  wakes);  or  they 
may  be  temporally  induced  by  the  presence  of  instabilities  or  gasdynamic 
discontinuities.  Such  instabilities  can  lead  to  large  scale  coherent 
structures  or  small  scale  structures  we  may  more  properly  term 
"turbulence"— either  electromagnetic  and/or  fluid  mechanical  (electrothermal, 
magneto -acoustic,  hydromagnetic  instabilities;  arcs). 

In  the  present  work  we  examine,  some  of  the  detailed  structure 
of  magneto  gasdynamic  flows  encountered  in  pulsed  magnetohydrodynamic 
(MHD)  generators  [l],  [2],  [3].  We  do  so  with  an  eye  to  revealing  the 
nonuniformities  which  may  be  significant  in  the  operation  of  actual  devices. 

STD  Research  personnel  have  devoted  considerable  research  to  the 
detailed  prediction  of  the  phenomena  in  and  performance  of  combustion- 
driven  and  closed-cycle  nonequilibrium  MHD  generators.  An  important 
aspect  of  the  detailed  predictive  approach  is  that  important  effects  were 
predicted  and  exhibited  before  they  were  observed  in  actual  experiments; 
hence  they  serve  as  an  important  guide  to  the  fundamental  design  of  such 
machines.  These  effects  include  inter -electrode  shorting  and  breakdown 
[4],  Lorentz  force  generated  velocity  profile  overshoots  (and  correspondingly 
enhanced  wall  friction  and  heat  transfer)  [5],  tranverse  electrode  arcing 
[6],  and  MHD  driven  secondary  flows  [7].  In  the  present  work  we 
describe  some  initial  efforts  to  bring  this  same  predictive  approach  to  bear 
upon  pulsed  MHD  generator  flows.  Because  of  a prime  requirement  to 
minimize  computation  costs,  these  studies  were  conducted  with  reduced  and 
simplified  versions  of  the  full  STD  Research  codes  in  which  full  multi- 
dimensional effects,  kinetic  and  diffusional  effects,  and  detailed  turbulence 
modeling  are  included.  Nevertheless  they  reveal  significant  effects  which 
the  more  detailed  code  implementations  can  subsequently  elaborate. 
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Dense  plasmas  which  serve  as  the  working  fluid  for  a class  of  these 
devices  pose  important  questions  as  to  the  precise  value  of  the  equilibrium 
electrical  conductivity  which  exists  at  attainable  explosion  generated  states 
(plasma  pressure  10-100  k bar,  plasma  energy  5-10  eV),  Various  theoret- 
ical estimates  show  that  the  highly  compressed  plasma  conductivities  may  1 

differ  by  nearly  a factor  of  2 from  the  classical  Spitzer  Harm  values  [5]. 

As  a contrast  to  these  dense  plasma  corrections  in  equilibrium 
properties  we  will  show  that  in  the  dynamical  or  flow  state  important  internal 

structure  exists  in  both  unsteady  shocked  MHD  flows  as  well  as  steady  state  > 

high  Magnetic  Reynolds  Number  channel  flow.  The  flow  nonuniformities 

associated  with  these  structures  induced  by  the  dynamical  interaction  of 

the  working  fluid  with  the  external  magnetic  field  can  alter  performance  by 

an  order  of  magnitude  from  that  for  uniform  flow. 

We  consider  two  illustrations  of  such  magnetogasdynamic  non- 
uniformities. The  first  is  for  conventional  explosion-driven  MHD  generators. 

Here  we  reveal  the  detailed  structure  of  the  deceleration  of  the  conducting 
slug  formed  by  the  compression  of  plasma  between  the  explosion-driven 
shock  wave  and  its  following  contact  surface.  Such  structure  includes 
electrothermal  and  magnetic  contraction  of  the  current  distribution  within 
the  slug  and  the  corresponding  nonuniform  deceleration  and  Joule  heating 

of  the  slug.  The  second  illustration  is  for  steady  high  Magnetic  Reynolds  I 

Number  channel  flow.  Here  we  exhibit  the  nonuniform  manner  in  which 
the  external  magnetic  field  is  diffused  into  the  plasma  and  the  corresponding 
nonuniform  deceleration  of  the  flow. 

2.  UNSTEADY  FLOW:  THE  SHOCK  TUBE  MHD  GENERATOR 

2 . 1 Characterization  of  Nonuniformities 

The  idealized  shock  tube  MHD  generator  [l]  is  shown  in  Fig. 

2-1.  An  appropriately  shaped  chemical  explosion  charge  is  detonated  into 
a gas  which  is  compressed  and  heated  behind  the  driven  shock  wave  and 
the  contact  surface  which  separates  the  fill  gas  in  the  tube  and  the  products 
of  explosion.  The  high  velocity  conducting  slug  so  formed  passes  through 
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DIRECT  EXPLOSION  DRIVEN  MHD  GENERATOR 
IDEALIZED  SYSTEM 


HIGH  TEMPERATURE,  HIGH  VELOCITY 


UNIFORM  PLASMA  FILLS 
ENTIRE  REGION  FROM  SHOCK 
TO  CONTACT  SURFACE 


8-2851 


Fig.  2-1.  Idealized  shock  tube  MHD  generator  flow 


3 


the  magnetic  field  generating  a uniform  current  over  its  breadth  [lj,  [9] 
and  is  uniformly  decelerated. 

In  Fig.  2-2,  we  exhibit  some  of  the  nonuniformities  which  will 
be  present  in  such  a flow.  The  shock  front  of  the  conducting  zone  will  not  be  a 
sharply  delimited  region,  but  will  in  fact,  be  dispersed  by  two  principal 
effects.  First,  kinetic -relaxation,  radiative  processes,  and  diffusionally 
induced  electric  fields  will  smear  the  shock  front  [10],  [ 1 1 ] . For 
conditions  typical  of  chemical  explosion  drives,  the  shock  front  kinetic 
dispersion  is  of  the  order  of  0.1-  1.0  cm.  [l0],[ll].  It  is  important  to 
note,  however,  that  careful  consideration  of  all  diffusional  processes  in 
a comprehensive  Ohm's  law[l3]  indicates  significant  effects  which  have 
not  previously  been  considered.  For  example.  Hall  effects  are  normally 
neglected  in  such  flows  because  of  the  high  pressures  involved;  however, 
it  is  possible  to  show  that  the  current  induced  by  the  diffusional  term 
Vpe  X B in  the  Ohm's  law  [l 3 ] is  in  the  same  direction  as  the  generator 
current  and  for  very  strong  shocks  has  the  limiting  value 

w 

J(Vpe  «)'*t 

where  e is  the  electronic  charge,  ug  is  the  electron  cyclotron  frequency, 

and  Q . is  the  electron-ion  cross  section.  Such  a current  is  of  the  order 
_ei 

of  10^  kA/m^  and  hence  makes  an  important  modification  to  the  shock  front 
of  a hydromagnetic  flow.  The  second  dispersal  mechanism  for  the  slug 
front  is  the  Lorentz  force  itself. 

A second  major  nonuniformity  in  such  a hydromagnetic  flow  is 
the  nonuniformity  induced  in  the  current  distribution  by  electrothermal 
heating  and  magnetic  contraction.  Since  shock-generated  conducting  slugs 
in  generators  of  length  0.1-1  m and  applied  magnetic  fields  of  1-10  T can 
be  significantly  Joule  heated,  the  uniform  discharge  hypothesized  for  the 
ideal  case  of  Fig.  2-1  is  electrothermally  and  magneto -acoustically  unstable 
as  shown  by  Demetriades,  et  al.[l4]and  Oliver  [ 1 5 ] . Hence,  this  electro- 
thermally unstable  discharge  will  contract  either  into  one  or  several  current 
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Fig,  2-2.  Shock-tube  MHD  generator  flow  with  fluid  and 
electrical  nonuniformities 
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concentrations  or  large  scale  'arcs'  depending  upon  the  nature  of  initial 
nonuniformities  present.  We  emphasize  that  electrothermal  contraction  is 
a phenomenon  which  can  exist  in  low  as  well  as  high  magnetic  Reynolds 
number  flows.  In  the  case  of  high  Rm  flows  however,  another  contraction 
mechanism  is  available  which  preferentially  drives  the  current  distribution 
to  the  front  of  the  conducting  slug.  This  is  the  hydromagnetic  contraction 
effect  which  concentrates  the  current  in  a progressively  thinner  'skin'  at  the 
slug  front.  Electrothermal  contraction  and  magnetic  contraction  can  of  course, 
operate  in  concert.  The  net  result  is  that  in  hydromagnetic  flow  in  contrast  to 
purely  gasdynamic  shock  driven  flow,  the  effective  length  of  the  conducting  siug 
is  not  the  region  between  the  shock  and  the  contact  surface,  but  rather  a 
hydromagnetic  structure  which  has  its  own  evolving  internal  nonuniformity. 

A third  class  of  nonuniformities  are  these  associated  with 
instability  and  mixing  at  the  contact  surface.  Such  anomolous  mixing  has 
been  observed  in  both  purely  gasdynamic  [ 1 6 ] and  hydromagnetic  flows. 

In  the  latter,  the  simplest  observation  is  that  alkali  seed  added  to  the 
explosive  significantly  affects  generator  performance.  In  the  idealized 
generator  of  Fig. 2-1  such  material  cannot  exist  in  the  conducting  zone. 

Lastly,  we  mention  wall  effects  which  include  boundary  layers 
developing  behind  the  shock  (are  they  truly  laminar  in  the  hydromagnetic 
case  as  evidence  seems  to  suggest  for  purely  gasdynamic  shocked  flow 
[17]?)  and  local  small  scale  arc-mode  transport  at  the  electrode  surfaces. 

2.  2 Mathematical  Description 

We  consider  a quasi-one-dimensional  description  of  a fluid 
moving  in  the  spatial  coordinate  x and  unsteady  in  time  t.  Average 
properties  (over  the  cross  section)  of  the  duct  describing  the  flow  are  its 
density  p,  internal  energy  e,  velocity  u,  pressure  p,  and  species  mass 
compositional  fraction  cft(o  = 1,2,  •••n).  From  these  we  select  the 
total  mass,  momentum,  energy  densities  [p,m  = pu,  e = p(e  +u2/2)],  as 
the  state  U (x,t)  specifying  the  flow  at  any  point  in  space  and  time: 


U(x,t)  = 


P 

m 

e 


(2.1) 


The  state  U.  is  augmented  with  the  compositional  vector  c . The  fluid 

a 

conservation  laws  are 

9tU  = -9xF  + 9x(e9xU)  + S - D (2.2a) 

9.c  + u9  c = 9 (e  9 c ) + c (2  2b) 

t a x a x'  x a'  a 

In  the  above,  £ represents  the  convect^d  fluxes  of  mass,  momentum,  and 
energy,  S and  D are  the  Lorentz  and  wall  loss  vectors,  e represents  a 
diffusion  coefficient  characterizing  the  mixing  processes  in  the  region  of 
the  contact  surface,  and  c is  the  chemical  source  term  for  species  a. 

In  terms  of  the  current  density  T,  electrict  field  E,  and  magnetic  field  if, 
the  Lorentz  vector  is 


(2.3) 


The  loss  vector  D has  the  form 


(2.4) 


, d&*  . 

where  is  the  rate  of  change  of  displacement  thickness  of  the  boundary 

layer  witl  respect  to  shock  coordinates,  and  t , qw  are  the  wall  shear 
stress  and  heat  flux,  respectively. 
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In  the  illustration  study  we  assume  that  the  kinetic  effects  are 
confined  to  the  relaxation  layer  at  the  shock  [we  do  not  consider  here  such 
phenomena  as  photoionization  and  gas  breakdown  in  regions  adjacent  to  the 
shock  [10] , [ 11  ] , [12]  which  can  occur  under  some  conditions.  ] Further,  for 
the  illustrations  to  be  presented,  we  assume  this  relaxation  layer  is  thin 
compared  to  the  overall  thickness  of  the  conducting  zone. 

The  electromagnetic  effects  are  described  by  the  Maxwell 

equations 


VXB  = p0J  8tB  = -V  X E (2.5) 

and  the  Ohm's  law  which  outside  the  region  of  intense  gradients  assumes 
the  form 


J=cr(E+uXB) 


(2.6) 


with  cr  = c(p,  T , c^) . The  forms  of  Eq.  (2.5)  allow  representation  of  E 
and  if  in  terms  of  scalar  potential  $ and  vector  potential  a as 


E = -V$  - 3tA 

B = V X T 


(2.7) 


For  the  simple  geometry  of  magnetic  field  in  the  z direction  and  continuous 
electrodes  perpendicular  to  the  y axis,  the  electrical  equations  reduce  to 


(V- 


E = - 

y 


V 


a.  a 

t y 


b = 9 A 
z x y 


(V-V  ) 

9.  A +u9A  = q 9 9a  - - — J.- 
ty  xy  xxy  n 


(2.8) 
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I 

V 


In  the  above,  V is  a uniform  voltage  over  the  range  of  the  continuous 
electrodes  and  Vj  is  the  voltage  drop  through  boundary  layers. 


Since  the  characteristic  electric  field  in  the  MHD  flow  is  of  the 

order  u X B we  may  estimate  the  importance  of  the  inductive  term  9ti* 

relative  to  u X B.  The  term  d A is  of  the  order  B.u  where  B is  the 

1 i 1 

field  induced  by  the  plasma  currents.  From  the  last  of  Eqs.  (2.8)  we 

find  B^~  (1-K)RmBz  and  V $ ~ KuBz  where  = p0<rua  is  the  Magnetic 

Reynolds  Number  and  a is  the  characteristic  length  of  the  conducting  zone. 

Hence  we  find  |3  A*|  is  of  order  (l-K)R  | u X B | . Thus,  self-induction 
l m ' 

effects  in  the  plasma  must  be  included  for  (l-K)R  > 1 

m 


The  electrical  equations  (2.8)  require  closure  with  an  external 
loading  condition.  If  we  integrate  the  last  of  Eqs.  (2.8)  in  x over  the 
electrode  extent,  we  obtain  the  total  current  as 


I = -A 


(O 


(V-Vd) 


( cr  u B ) 
' z' 


- A 


<°Vy> 


(2.9) 


where  A is  the  electrode  area  and  < ) indicates  an  average  over  the 
length  of  the  electrodes.  Defining  the  internal  resistance  of  the  generator 
as  R.  = h/(A ( cr ))  and  the  Faraday  voltage  Vp  as  VF=h(auB  >/<o->  we 
may  express  Eq.  (2.9)  as 


crd  A > 
t y' 


+ IR. 


■(V-Vd)  -vF 


The  external  circuit  equation  is 


L 

e 


RLI 


V 


(2.10) 

(2.11) 


where  Rl  and  Le  are  the  load  resistance  and  external  circuit  inductance 
respectively.  Eqs.  (2.10)  and  (2.11)  govern  the  evolution  of  the  global 
generator  current  I and  voltage  V. 
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Flow  Structure  Prediction  in  a Shock  Tube  Generator 

We  now  turn  to  an  application  of  the  foregoing  mathematical 
description  to  a typical  shock  tube  MHD  generator.  We  consider  a 
generator  of  approximately  1 m length  and  2.5  cm  cross  section  dimension 
whose  explosive  charge  is  an  RDX  composition  detonating  into  residual 
explosion  products  [2].  The  initial  and  operating  conditions  for  this 
machine  are  given  in  Table  2-1.  The  initial  condition  for  the  simulation 
is  taken  as  that  of  fully  detonated  explosion  products  at  detonation  pressure 
and  temperature  occupying  a volume  negligibly  different  from  that  of  their 
nondetonated  state  (Region  4)  and  residual  explosion  products  at  300°K  and 
10  Torr  pressure  (Region  1). 

In  the  interest  of  minimizing  costs  in  this  initial  study, 
simplified  thermodynamic  and  electrical  conductivity  functions  were 
utilized  rather  than  the  comprehensive  real  gas  thermochemical  and 
electrical  conductivity  algorithm  in  the  STD  Research  codes.  These  con- 
sisted of  a parameterized  specific  heat  ratio  y which  reflected  the  reactive 
nature  of  the  residual  explosion  products  (y  ~ 1.17)  and  simple  conductivity 
functions  allowing  for  a singly  ionized  species  interacting  with  a neutral 
background  (three  component  plasma).  Since  significant  multiple  ionization 
is  likely  in  these  plasmas,  the  simplified  functions  are  probably  not 
accurate  above  about  20,000°K.  Above  20,000°K,  this  simplified  con- 
ductivity is  suppressed  by  the  upper  limit  on  electron  number  density 
imposed  by  the  restriction  to  a singly  ionized  species.  This  conductivity 
function  will  underlie  the  structure  of  some  of  the  results  we  will  now 
describe . 

Since  the  computational  simulation  may  be  executed  with  or 
without  the  effect  of  induced  magnetic  field,  we  will  examine  the  differences 
between  flows  which  are  allowed  to  generate  an  induced  magnetic  field  and 
tnose  which  are  not.  In  the  first  sequence  of  results  in  Figs.  2-3a  through 
2-3d  and  Figs.  2-4a  through  2-4c,  the  flow  is  not  permitted  to  experience 
the  effects  of  an  induced  magnetic  field.  In  Figs.  2-3  we  observe  the 
evolution  of  the  flow  in  the  x-t  plane  with  the  amplitude  of  the  various 
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TABLE  2-1 


Nominal  Conditions  for  Shock  Tube 
MHD  Generator  (Ref.  [2  ]) 


Explosion  products  in  zone  4.  at 
pressure  p4,  temperature  T4  at  time  t = 0. 
Fill  gas  (residual  products)  in  zone  1 at 
pressure  pj,  temperature  Tj 


P4/P1 

2.2  X 107 

VT1 

17 

Ushock  ideal 

8236  m/s 

^shock  ideal 

27 

Bo 

1 Tesla 

Btube 

2. 54  cm 

^electrodes 

86  cm 

Load  Impedance 

4 X 10"3  ohm 

2%  potassium  seed  in  fill  gas 


8-2853 
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Evolution  of  current  distribution  in  shock  tube  MHD 
generator  for  conditions  of  Table  2-1.  Induced 
magnetic  field  suppressed.  Note  growing  electrothermal 
instability  in  slug  consisting  of  current  concentration  at 


LOG  ( PRESSURE ) g 5 g 5 g 2 

X Z X Z X 2 


F = L0GI PRESSURE) 

X = AXIAL  DISTANCE.  M 
T = TIME.  SEC. 

= 0. 

= 8-1 IE— 01 

= 0. 

= 1.48E-04 
= 2.43E+00 
= 9.78E+00 


8-2855 


Fig.  2-3b.  Pressure  evolution  for  conditions  of  Fig.  2-3a.  Note 
Increased  pressure  inside  conducting  slug  due  to 
decelerating  Lorentz  force 
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F = VELOCITY.  M/SEC 
X = AXIAL  DISTANCE.  M 
T = TIME.  SEC. 

X MIN  = 0. 

X MAX  = 1 .15E+00 
T MIN  = 0. 

T MAX  = 1 .50E-04 
F MIN  =-l .89E+01 
F MAX  = 7.03E+03 
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Fig.  2-3c.  Velocity  evolution  for  conditions  of  Fig.  2-3a.  Note 
maximum  deceleration  in  center  of  slug  and  velocity 
"overshoots”  at  slug  edges 


o 


• j O <D  CNi  ^ 

uj  a:  o o o o o 

C _» 

^ ^ i ii  ii  m ii  ii  ii 
Hdl-JXZXZX 

ii  ii  ii  3c  E 3c  £ S:  5 

l»-  X >-  X X t-  »-  ti_U- 


bc 

3 


ro  be 

'za 

• -M 

bo  <d  a 

•r*  Z.  rj 

t.2-o 

»M  3 0 

o 8 « 

3 « S 


So 

£ 60 

°a 

S rt 

n 2 

Rj  V 


.3  £ *3 

3 ft  ° 

g £ 

► a»  o 

« -m  *> 

If'0 

«i  o a 
h h O 
V bOT? 

Sft  H 

v 

-H  o 

« JO  > 

HZ  S 


CO 

I 

rg 


bO 

t4 

&H 


15 


Fig.  2 -4c.  Detailed  distribution  of  fluid  and  electrical  variables 
through  the  high  velocity,  high  temperature  region 
for  the  evolution  In  Figs.  2-3,  t = 140  ps.  Note 
"double-humped"  electrothermal  Instability  due  to 
single -Ionization  conductivity  law.  Full  real  gas 
effects  would  likely  lead  to  only  a "one -humped" 
instability. 
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fluid  and  electrical  variables  displayed.  In  Figs.  2-4  we  observe  the 
same  situation,  but  in  expanded  scale  so  as  to  observe  the  flow  structure 
present.  At  about  20  ps  the  conducting  slug  enters  the  generator  and 
produces  the  current  pulse  which  propagates  through  and  out  of  the  gener- 
ator by  about  160  ps.  The  overall  length  of  the  slug  grows,  but  the 
current  distribution  within  the  slug  has  its  own  evolving  internal  structure. 
Because  of  the  decelerating  Lorentz  forces,  the  pressure  within  the  slug 
is  actually  increased  over  that  for  a non-MHD  flow  and  an  expansion  wave 
follows  the  slug  as  shown  in  Fig.  2-3b.  Similarly,  the  velocity  within  the 
slug  is  decelerated  below  that  at  the  leading  and  trailing  edges  of  the  slug 
leading  to  velocity  "overshoots"  at  the  leading  and  trailing  edges.  (Fig. 
2-3c).  In  Fig.  2-3d  it  can  be  seen  that  there  is  a significant  growth  in 
the  average  temperature  levels  within  the  slug  due  to  the  Joule  heating 
present  (the  reader  should  bear  in  mind  that  the  actual  temperature  levels 
must  be  considered  approximate  due  to  the  use  of  simplified  thermodynamic 
functions).  In  Fig.  2-3a  and  Fig.  2-4a  it  can  be  seen  that  during  the 
course  of  its  evolution  a growing  nonuniform  current  concentration  is 
present.  This  is  because  in  the  range  of  temperatures  below  20,000°K 
(for  the  simplified  conductivity  model)  the  plasma  is  electrothermally 
unstable  [ 14 ] , [ 15 ] and  the  initial  nonuniformity  present  in  the  slug  shape 
is  growing.  Once  the  temperatures  exceed  20,000CK  in  the  center  of  the 
slug,  the  central  region  actually  becomes  electrothermally  stable  (because 
the  simplified  conductivity  functions  have  the  property  9<r/9T  < 0 in  this 
region).  Since  only  the  edges  of  the  slug  are  unstable,  these  edge  regions 
grow  and  we  see  a "double -humped"  growing  instability.  In  Figs.  2-5  and 
2-6  we  see  the  same  pulse  evolution  but  for  higher  interaction  corresponding 
to  Bq  = 3 Tesla.  Here,  the  electrothermally  unstable  slug  has  evolved 
greater  nonuniformities. 

For  real  gases  in  which  9<r/9T  is  a principally  positive  function  over 
the  entire  range  of  temperatures  and  pressures  encountered,  the  electro- 
thermal instability  will  likely  evolve  as  a single  growing,  contracting, 
intensifying  nonuniformity. 
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Fig.  2- 5b.  Velocity  evolution  for  BQ  = 3 Tesla  case  In  Fig.  2-5a 


In  Figs.  2-7  and  2-8  we  see  the  same  ease  shown  in  Figs.  2-3 

and  2-4  but  with  full  inclusion  of  induced  magnetic  field  effects.  The 

interaction  is  somewhat  higher  in  this  case  and  because  of  the  hydro - 

magnetic  contraction  of  the  current,  the  rear  instability  has  been  diminished 

and  the  leading  edge  instability  has  been  enhanced  and  driven  to  the  front 

of  the  slug  (Figs.  2-7a,  2-8a).  The  signature  of  the  induced  field  can  be 

seen  in  Fig.  2-7b  where  the  conducting  slug  has  evolved  a 1.48  Tesla  field 

ahead  and  an  0.52  Tesla  field  behind.  For  these  conditions  R =0.5. 

m 

In  Fig.  2-9  we  see  the  evolution  of  the  total  current  and  total 

~\  _ 3 

energy  produced  in  the  run.  Notice  that  at  this  4 X 10  12  loading  condition 

the  maximum  current  is  reached  near  the  end  of  the  run.  This  is  because 
the  current  density  levels  in  the  plasma  rise  with  time  due  to  Joule  heating 
and  axial  growth  of  the  slug.  The  4X10  12  loading  is  perhaps  closer  to 

short  circuit  rather  than  open  circuit  for  these  conditions  which  correspond 
to  those  of  [9].  Closer  to  open  circuit,  one  would  expect  less  Joule 
heating  and  the  maximum  current  point  would  occur  midway  in  the  pulse. 
Closer  to  short  circuit  the  maximum  current  would  occur  near  the  end  of 
the  pulse.  In  Figs.  2-10  we  display  some  of  the  experimental  results  for 
these  conditions  from  [9].  It  can  be  seen  that  the  shift  of  the  current  maximum 
to  the  late-time  portion  of  the  run  seems  to  correlate  with  the  degree  of  loading. 
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Fig,  2-7a.  Evolution  of  the  current  distribution  for  the  case  shown  in  Figs,  2-3  and 
204  but  with  induced  magnetic  field  effects  Included.  Note  how  inductive 
effects  have  concentrated  the  electrothermal  instability  at  the  front  of  the 
conducting  slug.  (See  also  Fig.  2-8) 
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Fig.  2-7b.  Evolution  of  the  magnetic  field  for  the  case 
1 of  Fig.  2-7a 


Fig.  2-9.  Evolution  of  the  total  current  and  energy  during 
the  run  shown  in  Figs.  2-3  and  2-4. 
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Fig.  2-10.  Experimentally  observed  total  current  pulse  shape 
from  the  experiments  of  Bangerter,  et  al,(  [5], 
Shift  of  maximum  current  to  late -time  portion  of 
pulse  appears  to  correlate  with  degree  of  loading. 


3.0  STEADY  FLOW;  MHD  CHANNEL  FLOW  AT  HIGH  MAGNETIC 
REYNOLDS  NUMBER 

3.1  Flow  Structure  Prediction  in  Hydromagnetic  Channel  Flow 

In  short  burst  MHD  generator  channel  flow  (timescale  100  ^s  ~ 

1000  ^s)  utilizing  explosively  generated  and  compressed  plasma  as  a stag- 
nation source  [3],  the  flow  is  very  nearly  a steady  channel  flow.  At  the 
plasma  conditions  of  electrical  conductivity,  length  scale,  and  velocities 
achievable  in  such  flows  both  the  interaction  parameter  I (ratio  of  Lorentz 
force  to  change  of  momentum)  and  magnetic  Reynolds  number  R^  are  large. 
For  the  hydromagnetic  flows  typical  of  the  generators  characterized  by  the 
operating  conditions  of  Table  3-1  we  find  the  magnetic  Reynolds  number  Rm 
to  be  about  20  and  the  interaction  parameter  I to  be  about  15/meter.  These 
conditions  indicate  that  the  deceleration  time  for  the  fluid  is  of  the  order  of 
the  diffusion  time  for  the  external  magnetic  field  (consisting  of  applied  and 
possible  self-excited  combinations)  to  penetrate  through  the  plasma.  One 
may  therefore  expect  very  significant  interaction  between  the  velocity  field 
of  the  flow  and  the  nonuniform  Lorentz  forces  associated  with  the  nonuniform 
magnetic  field  diffusion. 

In  Fig.  3.1  we  indicate  schematically  the  interaction  of  the  mag- 
netic field  with  the  high  Rm  channel  flow.  The  frontal  structure  (which  is 
not  significant  over  the  course  of  pulses  — 10  jis  and  which  we  will  not 
consider  further)  consists  of  axial  compression  and  bending  of  the  magnetic 
field  lines. 

Away  from  the  front  and  throughout  the  channel,  the  lines  are 
more  sharply  bent  through  the  magnetic  boundary  layers  and  resemble  the 
cross  plane  structure  shown  in  Fig.  3.1.  At  high  magnetic  Reynolds 
number  the  magnetic  boundary  layers  are  thin  and  transverse  gradients  in 
y,z  are  much  stronger  than  axial  gradients  in  x.  Correspondingly,  the 
generator  current  is  concentrated  in  the  magnetic  boundary  layers.  At 
the  channel  inlet  the  magnetic  field  is  most  strongly  excluded  from  the 
flow;  at  the  exit  the  field  has  penetrated  most  deeply.  We  now  ask  how 
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TABLE  3-1 

Conditions  for  Illustrative  High  Magnetic 
Reynolds  Number  Channel  Flow 


Generator  Inlet  Conditions 

<u> 

T 

P 

<r 


7 X 103  m/s 
2.8  X i04#K 
3 k.bar 

5 X 104  mho/m 


External  Magnetic  Field  (Bq) 
Generator  Cross  Section  (a) 
R 

m 

I 


10  T 
5 cm 
22 

15  m 


8-2869 
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Axial  and  Frontal  Structure  of  Magnetic  Field 
in  Hydromagnetic  Channel  Flow 


Cross  Plane  Structure  of  Magnetic  Field 
in  Hydromagnetic  Channel  Flow 

Fig.  3-i.  Magnetic  field  nonuniformlties  in  high  Magnetic 
Reynolds  Number  channel  flow 
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a hydromagnetic  flow  will  respond  in  its  velocity,  current,  Lorentz  force, 
and  power  generation  distributions  to  this  fundamental  nonuniformity 
situation  in  high  Magnetic  Reynolds  Number,  high  interaction  parameter  flow. 

3.  2 A Reduced  Mathematical  Model 


We  center  attention  on  the  axial  velocity  ux(x,y,z)  and  neglect 
the  effects  of  the  cross  flow  or  secondary  flow  velocity  field.  Since  I»  1, 
the  flow  is  dominated  by  electromagnetic  rather  than  fluid  mechanical 
viscous  effects. 

Under  conditions  of  I » 1,  R^  >>  1,  it  is  possible  to  show 

that  the  transverse  Lorentz  force  is  of  order  l/R  times  the  axial  Lorentz 

m 

force.  Hence,  cross  plane  momentum  is  weakly  generated  by  Lorentz 
forces  compared  to  axial  generation. 

The  primary  source  of  cross  plane  velocity  nonuniformities  is 
the  redistribution  of  axial  mass  flow  due  to  the  nonuniform  deceleration; 
this  is  not  an  insignificant  effect,  but  we  expect  a quantitative  rather  than 
qualitative  modification  to  simulations  in  which  this  effect  is  neglected. 

Consistent  with  our  neglect  of  cross  plane  momentum,  we  approxi- 
mate the  axial  pressure  gradient  3p/3x  as  3(p)/3x  where  ( ) indicates  an 
average  over  the  cross  section.  Mass  and  momentum  conservation  lead  to 

3x(pu>  =0  (3.1) 

(7xB) 

u3  u = - — 3 (p)  + 3t  (v3t  u)  + — (3.2) 

where  v is  the  viscous  diffusivity,  and  ^ = (y,z)  is  the  cross-plane  coordinate. 
In  the  case  of  laminar  flow,  v is  the  laminar  viscosity;  in  turbulent  flow 
v = aqA  where  q is  the  RMS  turbulence  velocity,  a is  a universal  structure 
parameter,  and  A is  the  turbulent  scale  function.  In  fact,  for  the  conditions 
of  the  illustration,  viscous  effects  represent  a very  small  effect  for  the 
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momentum  balance.  The  pressure  gradient  is  determined  by  averaging 
the  momentum  equation  with  (u  9u/9x)  * 0: 


9X<P>  = <UXB)X>  - A'1  £ rwd£ 


(3.3) 


where  A is  the  cross  sectional  area,  t is  wall  shear  stress,  and  Q ds 

W 9 J 

represents  an  integral  over  the  periphery  of  the  cross  section. 

Let  us  now  consider  the  electrical  description.  With  u = (u  ,0,0), 
Eqs.  (2.7)  yield 


ux  9XB  = -B  9xux  + x ( B ' V)  Ux-V  X (q  V X B" ) (3.4) 

where  x is  the  unit  vector  in  the  x direction  and  rj  = (p^o-)"1  is  the 

magnetic  diffusivity.  The  boundary  conditions  for  B are  as  follows.  As 

can  be  seen  from  Fig.  3-1,  specifications  of  the  boundary  data  for  IT  on 

the  channel  perimeter  requires  solution  of  the  external  (outside  the  channel) 

as  well  as  internal  problem  for  B . As  an  alternative  to  solution  of  the 

external  problem,  we  approximate  the  boundary  conditions  for  B B as 

v’  z 


W 

°n  z = ± J 


¥ (Bq  + AB)  exp(-\x)  y > 0 
± (Bq  + AB)  exp(-\x)  y < 0 


Bz=  B0  t 1 " exp(-X-x)  ] 


(3.5) 


w 

on  y = ± 7 


= Bn  + AB  exp(-Xx) 


(3.5) 
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i \ 47T  T)  Z 

where  \ j is  the  decay  term  for  the  fundamental  diffusion  mode 

in  a square  duct,  Bq  is  the  far-field  external  value  of  the  magnetic  field 
(totally  in  the  z direction),  and  AB  is  the  maximum  compression  of  the 
magnetic  field  at  the  generator  inlet.  It  can  be  seen  that  the  approximate 

model  (3.5)  possesses  the  correct  limiting  values  at  x — 1 - oo.  For  R » 1 

m ’ 

the  values  at  x — *■  0 are  also  close  to  the  full -field  exclusion  solutions. 

The  model  is  least  accurate  for  conditions  of  R =1 

m 

For  B the  boundary  conditions  are  B = ±u..~onz=±^r 
x x 2 2 

where  I is  the  total  current  per  unit  depth  flowing  in  the  generator  in 
the  transverse  direction. 

3. 3 Magnetic  Diffusion  and  Flow  Interaction 

We  now  consider  the  predictions  for  flow  in  the  generator  channel 
of  Fig.  3-1  under  the  conditions  of  Table  3-1.  For  these  conditions  we 
set  the  generator  current  at  the  nominal  value 

I = (1  -K)  w ( <r  u B ) 

with  K = l/2  and  w = the  channel  width.  We  assume  that  no  field  has 
penetrated  the  plasma  at  the  electrode  edge  x = 0.  (In  general,  one  will 
expect  a starting  magnetic  boundary  layer  thickness  at  the  channel  entrance 
dependent  upon  the  plasma  nozzle  geometry  and  the  fringe  field  location.) 

The  distribution  of  the  transverse  current  and  the  magnetic 
field  component  Bz  for  various  axial  stations  in  the  channel  are  shown  in 
Fig.  3.2.  (There  is  an  axially  directed  eddy  current  Jx  associated  with 
the  fields  Bz  and  B^  which  is  not  shown  here.)  At  x = 1.4  cm  the  mag- 
netic boundary  layers  are  exceedingly  thin  3 mm).  They  have  grown 
to  1.4  cm  at  x = 14  cm  and  2.9  cm  at  x = 28  cm.  The  inward  diffusing 
magnetic  field  component  Bz  is  initially  concentrated  along  the  walls 
parallel  to  the  external  magnetic  field  (electrode  walls)  (see  Fig.  3- lb), 
but  the  z component  on  the  side  walls  builds  up  as  the  flow  proceeds  down 
the  channel.  The  transverse  current  is  initially  concentrated  along  the 
sidewalls  and  diffuses  inward  from  these  walls. 


34 


Y 


In  Fig.  3.3  we  observe  the  impact  of  these  distributions  of 
current  and  magnetic  field  on  Lorentz  force  and  axial  velocity.  In  the 
early  stages  of  the  flow  evolution  the  maximum  axial  Lorentz  force  is  in 
the  corners.  This  is  because  the  maximum  field  exists  on  the  electrode 
walls  parallel  to  the  z axis  and  maximum  current  J exists  on  the  sidewalls 
parallel  to  the  y axis.  The  maximum  product  of  the  two  therefore  exists  in 
the  corners.  Throughout  the  entire  evolution,  the  maximum  Lorentz  force 
is  concentrated  in  the  magnetic  boundary  layers  along  the  sidewalls. 

The  impact  on  the  velocity  distribution  is  strong  and  significant. 

The  initially  flat  velocity  profile  is  rapidly  decelerated  in  the  sidewall 
layers.  By  12  cm  the  flow  is  actually  reversed  in  the  corners  and  the 
regions  of  reverse  flow  spread.  Because  of  the  deceleration  of  the  wall 
regions  by  the  inward  diffusing  Lorentz  force  the  core  flow  is  accelerated. 
From  an  initial  value  of  7 km/ s the  core  flow  is  accelerated  to  nearly 
10  km/s  by  x = 28  cm.  Since  there  is  little  Lorentz  force  in  the  core 
region,  this  high  velocity  region  is  not  effectively  coupled  by  the  generator. 

It  should  be  noted  that  once  reverse  flow  occurs  (stalled  MHD 
boundary  layers)  the  present  prediction  of  the  actual  velocity  field  in  the 
wall  layer  region  cannot  be  considered  correct  since  upstream  feedback 
is  prohibited  by  the  method  of  calculation.  The  present  calculation  is  only 
definitive  in  its  prediction  that  stall  will  indeed  occur.  In  actuality  we 
conjecture  that  the  flow  under  these  conditions  will  appear  as  shown  in 
Fig.  3-4  in  which  recirculation  eddies  are  trapped  within  the  boundary 
layer  region  of  intense  Lorentz  force  while  the  core  flow  jets  through 
outside  the  region  of  the  eddies.  In  such  a stalled  flow,  there  is  little 
Lorentz  force  coupling  to  the  high  momentum  jet  in  the  central  region  of 
the  channel. 

As  a demonstration  that  stall  is  not  a peculiarity  of  the  corner 
region,  we  show  the  Lorentz  force  and  velocity  profiles  at  x = 21  cm  for 
a higher  interaction  situation  (cr  = 2 X 10  mhc/m)  but  otherwise  conditions 
identical  to  Table  3-1  (Fig.  3-5).  In  Fig.  3-5b  we  see  that  the  entire 
sidewall  has  reversed  flow  under  these  conditions. 
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Fig.  3-3.  Lorentz  force  and  axial  velocity  distributions  at 

three  axial  stations  in  high  Magnetic  Reynolds  Number 
channel  flow.  Wall  layers  are  sharply  decelerated 
with  wall  layer  flows  actually  reversing.  Core  flow 
is  accelerated. 


Mf3X  = 5.554E+10 


Fig. 
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3-5a. 


Lorentz  force  distribution  at  x = 21  cm 
case  (cr  = 2 X 10^  mho/m)  but  otherwise 


corresponding  to  higher  interaction 
conditions  as  in  Table  3-1 
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In  Fig.  3-6  we  show  the  distribution  of  power  production  -(  J • E ) 
over  the  cross  section.  Power  is  only  produced  in  the  mid-wall  layer 
regions  and  is  actually  consumed  in  the  corner  and  sidewall  layers  where 
reversed  flow  takes  place.  The  overall  power  extraction  per  unit  length, 
dP/dx,  is  defined  as  the  integral  of  the  power  density  over  the  channel 
cross  section: 

A 

The  ideal  power  extraction  for  uniform  flow  is 

(IS)  ■=  K(i-K)<r„^A 

ideal 

In  Fig.  3-7  we  exhibit  the  ratio  of  these  quantities  for  the  conditions  of 
Table  3-1.  It  can  be  seen  that  for  a channel  15  cm  long,  only  10%  of 
the  ideal  power  is  available  and  for  a 28  cm  channel  only  27%  is  avail- 
able. These  results  must  be  tempered  by  the  observation  that  beyond 
12  cm  the  boundary  layers  have  begun  to  stall;  it  is  not  likely  that  power 
extraction  will  be  enhanced  by  realistic  inclusion  of  the  recirculation  zones 
under  stalled  conditions.  Hence,  the  power  extraction  levels  beyond 
x = 12  cm  will  likely  not  be  reached. 

It  is  appropriate  to  mention  at  this  stage  that  previous  experi- 
ments at  Hercules,  Inc.  on  pulsed  MHD  generators  failed  to  correlate  with 
the  conductivity  computed  by  Hercules,  Inc.,  and  their  consultants.  The 
conductivity  computed  by  Hercules,  Inc.  was  too  high. 

However,  when  the  STD  Research  theory  and  conductivity  cal- 
culations were  used  with  the  same  inputs,  the  calculated  conductivity  was 
approximately  40%  lower  and  there  was  excellent  agreement  with  experi- 
ments . 
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x = 14  cm 


3-6.  Power  extraction  density  evolution  for  the  case  In 
Figs.  3-2  and  3-3.  Power  Is  produced  In  m laye 
and  core  regions  but  consumed  In  near  wall  regloi 


Fig.  3-7,  Ratio  of  power  extraction  per  unit  length  to  ideal 
power  extraction.  Wall  layers  stalled  over  most 
of  channel  length 
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